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We consider the growth of a polymer layer on a flat surface in a good sol- 
vent by in-situ polymerization. This is viewed as a modified form of diffusion- 
limited aggregation without branching. We predict theoretically the forma- 
tion of a pseudo-brush with density (j>(z) oc z -2 / 3 and characteristic height 
H oc i 3 . These results are found by combining a mean- field treatment of the 
diffusive growth (marginally valid in three dimensions) with a scaling theory 
(Flory exponent v = 3/5) of the growing polymers. We confirm their validity 
by Monte Carlo simulations. 

PACS. 82.35+t,68.70+w,82.65-i 



Polymer brushes (layers of linear polymer chains end-grafted to a surface) find use in 
the stabilization of colloidal particles in solution One way to form these layers is 

by adsorption of pre-existing polymers having a functionalised group on one end, but this 
process turns out to be extremely slow for high grafting densities Q. An alternative is the 
polymerization of a brush in situ; brushes formed this way have been reported, but not 
characterized in detail j|. In this paper we study an idealized model of in-situ growth of 
a layer, in which a flat impenetrable surface is initially densely covered with initiator sites. 
From these seeds, linear unbranched chains can grow by irreversible polymerization, with 
monomers added only at chain ends. The chains so formed are treated as flexible coils: they 
can relax their structure in response to excluded- volume interactions with their neighbours. 
This process sets a dynamical time scale; we restrict attention to slow enough growth, such 
that local equilibrium is maintained at all times. We also assume a constant diffusivity for 
the free monomers and a constant, infinitesimal flux of incoming random walkers. (For a 
finite, or a time-dependent, incident flux, our description breaks down, as discussed later.) 

The above-formulated problem of in situ polymerization is closely related to that of DLA 
(Diffusion-Limited Aggregation ||) without branching [H^|[9j. There is no satisfactory 
mean-field fomulation for ordinary DLA Hj,[L(J], but a continuum field approach does exist 
to describe the growth of a forest of needles by diffusive deposition, with adsorption on the 



needle-tips only || . It was found by simulation || that the mean-field calculation is rather 
accurate in d = 3 dimensions, and later realized || that this is the upper critical dimension 
for the problem, so that only logarithmic corrections to the mean-field results should apply. 
In what follows, we generalize the needle growth model to a much wider class of processes, 
in which there is an arbitrary power-law exponent v defined by <f)(z) ~ £{z) 1/v - d where 
4>(z) is the density at height z and the associated transverse separation between linear 
growing objects. (This is the usual relation, for an object of fractal dimension 1/z/, between 
the density and a "blob size" which is locally given by the separation £,(z).) For in-situ 
polymerization in a good solvent, v is simply the usual Flory exponent v = 3/5 |TIJ. In this 
case, we predict the formation of a novel power-law "pseudobrush" [13[] in which chains are 
extremely polydisperse but, nonetheless, strongly extended from the surface: the density of 
grafted material decays with a power law of height, (j)(z) oc 2~ 2 / 3 , and the spacing between 
chains at height z varies as £ ~ z 1 ^ 2 . We have tested these predictions by Monte Carlo 
simulation using the bond- fluctuation model (BFM) ]TJj. Good agreement is obtained in 
d = 3, which remains the upper critical dimension for any v (as may be confirmed by a 
simple reworking of Krug's argument M). 
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FIG. 1. Sketch of diffusive growth of a polymer layer by in situ polymerization. 



The continuum field description of the growth can be formulated quite simply. We define 
three coarse-grained densities (each expressed as a local volume fraction): (i) that of the 
aggregated monomers (p(z,t); (ii) that of active end monomers p(z,t); and (hi) that of free 
monomers u(z,t). These are coupled by three differential equations, as follows. Firstly, we 
suppose the polymerization is a local second-order process: 

d t (j) = kup (1) 

where k is a reaction constant. Secondly, conservation of the particles implies: 

d t( j) = -d t u + DAu = DAu (2) 

where D is the diffusivity of free monomers. The second equality embodies the "adiabatic" 
approximation, valid in the quasi-static limit of infinitesimal flux M. Thirdly we have 
a "constitutive relation" between local monomer density and blob size, or equivalently 
between and the end density, p: 

p = -ld z ^ e (3) 
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Here / is of order of the monomer size, b. This relation arises because the number of 
chains M{z) per unit area of surface surviving to height z is (by simple geometry) of order 
g^ z y-d ^(i-d)/(i/i/-d) ) an( j the end density is —dMjdz. Hence e obeys e = (d—l/v)/(d—l). 
Within this continuum approach, the exponent e contains all required information about 
local chain structure. Adopting Flory's estimate, v — 3/(2 + d), we have e = 2/3 (which is, 
by accident, independent of d). 

Next we consider the boundary conditions. Within the continuum model, as we show 
below, the density of adsorbing chain ends is (formally) divergent near the wall, so the 
chance of an incident particle actually reaching the wall is negligible. We thus require that 
the density of the free monomers, and its gradient, vanish there: 

u(0,t) = du(0,t)/dz = 0. (4) 

Secondly, the boundary condition on u at infinity is that appropriate to a (small) constant 
flux joo of incoming monomers (each of volume b ): 

Jim du{z, t)/dz = b d Joo /D = 9oo . (5) 

This ensures (in the adiabatic limit) that the mass of the layer grows linearly in time. For 
convenience, we may set 1/k = I = 1 by choosing these as the units of time and length. 
Two parameters then remain, the reduced flux g^ and the reduced diffusivity D/kl 2 = D. 
The latter may itself be set to unity by the rescalings = D e (p,p = Dp,u = D e ~ x u and 
goo = D e ~ l g~ 00 . In what follows we suppress the tildes (i.e. we choose D — 1) but restore 
the scale factors later as required. 



We now define a characteristic layer height by the first moment ETJ Hit) 



J zcf)dz/ J (pdz, and seek solutions of Eqs. ([!]-§) in terms of the scaled variable x = z/H(t): 

4> = z~ a f(x) ; p = z~°h(x) ; u = z~> g(x) (6) 

Consistency of this scaling ansatz with the governing equations and the boundary conditions 
fixes the exponents a, (3, 7 as we now show. Firstly, we eliminate d t <p from Eqs. ([I], ||]) to ob- 
tain up = d\u. This forces the end-density to decrease with an exponent (3 = 2. Secondly, the 
constant flux boundary conditions, eqns. (|4]||) impose an exponent 7 = 1 (with g{oo) = goo 
and g(0) = 0). Finally, we require on physical grounds that the function f(x) approaches 
a finite limit at the origin, and that /(oo) = 0. The constitutive relation eq. (|3D, between 
monomer and absorber density then reads ez~ l3 h(x) = z~( 1+a / e ' > /(x) 1 / 6 (a — a; (In f(x))') 
where the prime denotes a derivative. Equating powers of z we obtain f3 = 1 + a/e; within 
mean field theory (j3 = 2), this implies a = e which fixes our final exponent. It follows imme- 
diately that £ oc z l ^ d ~ l \ and hence £ oc \fz in 3 dimensions; this result is independent of e. 
In contrast, the total mass M = / °° cj)(z)dz obeys M = •r]H 1 ~ e where 77(e) = / °° x~ e f(x)dx. 
Under constant flux conditions, M = goot, and one has accordingly Hit) = (goot/i]) 1 ^ 1 ^^. 
This completes, at the scaling level, the generalization of the mean field theory for needle 
growth (e = 1) || to the case of arbitrary e. Note that, at this level, the only difference 
between the needle model and the case of general e lies in the fact that, when a particle 
is adsorbed, the adsorbing end moves forward by one monomer size for needle growth but 
less than this for the polymer case. Since only the active ends participate in screening, the 
scaling p ~ z~ 2 remains as for needles, and (3 = 1 + a/e then follows from Eq.3. 
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To test these predictions, we have performed detailed bond-fluctuation model simulations 
(for the method, see @]). In this athermal model, "local jumps" are made by choosing one 
monomer at random (regardless of whether or not it is part of a chain) and attempting 
to jump one lattice unit in a randomly chosen basis direction. The attempt is accepted 
if excluded volume and chain connectivity requirements are satisfied. If a free monomer 
happens to touch the wall at z = a (a is the lattice spacing) it adsorbs and becomes 
reactive. A free monomer can polymerize at a reactive chain end whenever the difference 
vector of the two monomers is an allowed intrachain bond. The simulation box is of size 
L x = L y = 100a, L z = 500a with periodic boundary conditions in x— and ^/-directions and 
impenetrable walls at z = a and z = L z . The flux used was = 1.3 ■ 10~ 7 /ra 2 where 
t is the time for one Monte Carlo step per monomer. No hydrodynamic interactions are 
taken into account. In the adiabatic limit of interest, chains are fully relaxed and the precise 
dynamics should be irrelevant. 
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FIG. 2. The density profiles for the grafted monomers <f>, the end monomers p and the free 
monomers u versus distance z from the wall. An incident flux of joo = 1.3 • 1(T 7 /™ 2 was imposed 
over a period of 4 • 10 6 r. We compare the densities with the power law predictions a = e = 2/3 
(solid line), (5 = 2 (dashed line) and 7 = 1 (dashed-dotted line). Near the wall the density of free 
monomers increases with distance as u(z) = zg(x) ~ z 3 ~ e (dotted line) corroborating 7=1 and 
the asymptotic power law solution of eq.(0) for small x jl5|| . The free monomer density (calculated 
from the incoming free monomer flux) is very small and is vertically shifted. 

The density profiles for the grafted monomers 0, the end monomers p and the free 
monomers u are shown in fig. and confirm broadly the scaling predictions. In fig. ^| the 
scaling functions g(x) and f(x) are shown. The direct verification of the scaling prediction 
for the monomer density is less successful than for the free monomer density. The accuracy 
of f(x) is roughly indicated by the size of the wiggles, which apparently arise from statistical 
noise: the accuracy decreases with both x and t, as the density becomes very low. Also, 
the interactions of the first blobs with the wall and fluctuations near the outer brush-edge, 
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omitted from our course-grained description, imply two additional length scales (the sizes 
of the first and the last blob of the layer) which, judging from fig. [3] (b), are non-negligible 
even at our longest times. 




FIG. 3. The scaling functions (a) g(x) and (b) f(x) for serveral times. The solid line indicates 
the numerical solution of eq. (0), with a discontinuity at x c = 12, obtained by a Runge-Kutta 
method, (a) The data collapse justifies the scaling ansatz eq. (|6|) with an exponent 7 = 1. (b) For 
t = 10 6 r the grafted chains are not stretched. As the layer becomes more brush-like the curves 
merge. 

We now return to the scaling functions f(x), g(x) and h(x). Defining a dimensionless 
parameter /i = the mean-field equations (|l|-|5|) can be reduced to the form: 

h(x) g(x) = -fix^f = g(x)f(x)^ (l - 

= x(2g'(x) + xg"(x)) (7) 

These equations (combined with the boundary conditions) determine in principle the scal- 
ing functions f,g, h; they represent a powerful extension, to arbitrary e < 1, of the existing 
results for needle growth (e = 1). For the latter case, an exact solution of the equations 
has in fact been obtained ||. For e < 1, however, we have found numerically (and shown 
analytically ||15|| ) that no acceptable solutions satisfy the boundary conditions on /. This 
suggests that either the scaling ansatz is wrong (which seems unlikely, in view of the simu- 
lation data shown above), or that new physics intervenes. The latter may certainly be true 
within the outermost blob of the polymer layer, in which both the mean-field averaging, and 
the constitutive relation (3), must break down. Therefore we postulate the existence of a 
step discontinuity in f(x) at some value x c , with / = for x > x c . (This sharp step, and 
corresponding features in h and g, would in reality be spread over the a scale ~ ff(t) 1 / 2 , 
the size of an outermost blob.) This ansatz is enough to determine the scaling functions 



f,g, h uniquely, including the value of x c p5| , although in practice it is hard to obtain an 



accurate numerical value; our best estimate is x c = 12 ± 2. Numerical results from eq. (|7|) 
for f(x) and g(x) are shown alongside the simulation data in figs. Since H(t) can be 
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extracted directly from the data, there is no ambiguity on the horizontal scale, but the 
parameter values (k,l...) are not exactly known for the simulation, so g(x) is normalized 
vertically by the plateau value at x = oo. On this basis, agreement is good for g(x) and 
reasonable for f(x). The vertical scaling could be avoided by estimating parameter values 



in the simulation; agreement remains acceptable [15] . 

We now discuss the validity of the adiabatic approximation {d t u <C DAu) used in eqn. (0) 
(see J7]). By substituting the scaling results derived above (restoring dimensional factors), 
we find that this assumption is selfconsistent for short times, but breaks down when the layer 
reaches a crossover height h* « (D / kl) e ^ 1+e ^ g o~^ 1+e ' ) - This can be made arbitrarily large 
by taking the limit of small incident flux and/or rapid diffusion. For experimental purposes, 
it may anyway be more realistic to impose, instead of a constant flux, an asymptotic density 
of the free monomers in the bulk. This would entail a second characteristic length (the 
depletion length £d ~ \J~Di of the incident diffusers). So long as h <C £d, the incident flux 



is joo oc Uoo^D/t, and the layer mass M increases as \ft rather than linearly. Expressed in 
terms of M, however, the structural scaling exponents a and (3 should be unaffected. Once 
again, the theory will break down at long times. 

We have shown how the case of diffusion-limited in-situ polymerization of a polymer 
brush can be viewed as one of a large class of problems, related to needle-growth models 
BllJi)!]) each characterised by an exponent e; e ~ 2/3 for polymerization in a good solvent 
whereas e = 1 recovers the needle growth limit. We found that the scaling exponents 
for the density, end density and free monomer density obey a = e, (3 = 2,^ = 1. These 
scaling predictions are consistent with detailed Monte-Carlo simulations based on the bond- 
fluctuation model. This is not obvious because of the coarse-graining implicit the continuum 
description - which appears to be especially delicate in the outermost blobs of the layer (cf. 
our discussion of f(x) near x = x c ). An important extension of our work is to the case of 
polymerization onto a spherical colloidal particle ||. We intend to discuss this in a future 
paper, along with results for flat layers in d = 2 ||15||. 
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